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ABSTRACT 

An analytic model of rotating hot spot is proposed to explain kilohertz quasi- 
periodic oscillations (kHz QPO) in X-ray binaries, which is based on the magnetic 
coupling (MC) of a rotating black hole (BH) with its surrounding accretion disc. The 
hot spot in the inner region of the disc is produced by energy transferred from a 
spinning BH with non-axisymmetric magnetic field. The frequency and energy band 
of the hot spot turns out to be related to six parameters, by which the strength 
and the azimuthal angular region of the bulging magnetic field on the BH horizon, 
the mass and spin of the BH, and the power-law index of the magnetic field varying 
with the radial coordinate on the disc are described. In addition, the correlation of 
the fluctuation of the bulging magnetic field with the widths of QPO frequency is 
discussed. 
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1 INTRODUCTION 

With the existence of a magnetic field connecting a Kerr 
black hole (BH) to its surrounding disc, energy and angu- 
lar momentum can be transferred from the BH to the disc 
(Blandford 1999; Li 2000; Li 2002a, hereafter Li02a; Wang, 
Xiao & Lei 2002, hereafter WXL), which can be regarded 
as one of the variants of the Blandford-Znajek (BZ) pro- 
cess proposed two decades ago (Blandford & Znajek 1977). 
Henceforth this energy mechanism is referred to as the mag- 
netic coupling (MC) process. Recently Wilms et al. (2001) 
found that XMM-Newton observation of the nearby bright 
Seyfert 1 galaxy MCG-6-30-15 reveals an extremely broad 
and redshifted Fe Ka line indicating its origin from the very 
most central region of the accretion disc, and they suggested 
that the rotating energy of a BH is extracted by the mag- 
netic field connecting the BH or plunging region to the disc. 
Later Li pointed out that the magnetic coupling of a spin- 
ning BH with a disc can produce a very steep emissivity 
with index a = 4.3 ~ 5.0, which is consistent with the above 
XMM-Newton observation (Li 2002b, hereafter Li02b). This 
result can be regarded as the observation signatures of the 
existence of the MC process. Very recently we obtained the 
same result in a more detailed MC model, which is consis- 
tent with the XMM-Newton observation with reference to a 
variety of parameters of the BH-disc system (Wang, Lei & 
Ma 2003, hereafter WLM). 

Another relevance of MC model to the observation is 



kilohertz quasi-periodic oscillations (kHz QPO) in X-ray bi- 
naries. As argued by van der Klis (2000), kHz QPOs in X- 
ray binaries probably originate from the inner edge of an 
accretion disc with a BH of stellar-mass, since millisecond is 
the natural timescale for accretion process in these regions. 
Recently Li suggested that the non-axisymmetric MC of a 
spinning BH of stellar-mass with a disc might be used to 
explain the kHz QPOs (Li 2001, hereafter LiOl). Motivated 
by Li's suggestion we propose an analytic model of a ro- 
tating hot spot to explain the kHz QPOs, where the hot 
spot arises from the energy transferred to the inner region 
of the relativistic disc by the closed field lines of a non- 
axisymmetric magnetic field on the BH horizon. Six param- 
eters are involved in our model. Among them three are used 
to describe the strength and azimuthal angular region of the 
bulging magnetic field on the BH horizon, another two are 
used to describe mass and spin of the Kerr BH, and one is 
the power-law index of the magnetic field varying with the 
radial coordinate on the disc. We give a detailed discussion 
on a rotating hot spot for fitting kHz QPOs in terms of these 
parameters. 

This paper is organized as follows. In Section 2 we 
show that the position of a hot spot is determined by non- 
axisymmetric magnetic field with a mapping relation be- 
tween the poloidal angular coordinate on the horizon and 
the radial coordinate on the disc. It turns out that the hot 
spot is limited very close to the inner edge of the disc for the 
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Figure 1. The poloidal magnetic field connecting a spinning BH 
with its surrounding disc 
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Figure 2. Azimuthal profile of the non-axisymmetric magnetic 
field on the BH horizon 



appropriate values of the BH spin and the power-law index 
of the magnetic field on the disc. In Section 3 we prove that 
the hot spot can reach the energy level of emitting X-ray, 
provided that the bulging magnetic field on the horizon is 
strong enough to reach ~ W s Gauss. Furthermore the rotat- 
ing frequency of a hot spot is expressed in terms of Keplerian 
angular velocity at the place where r 2 F is maximized, where 
F is the sum of radiation flux due to disc accretion and that 
arising from the MC between the spinning BH and the disc. 
In Section 4 we fit kHz QPO widths by assuming the fluctu- 
ation of the bulging magnetic field on the horizon. Finally, 
in Section 5, we give a summary for our model. 

Throughout this paper the geometric units G = c = 1 
are used. 



2 RADIAL POSITION OF ROTATING HOT 
SPOT 

In order to facilitate the discussion of the rotating hot spot 
we make the following assumptions. 

(i) The magnetic field is assumed to connect the BH 
with the surrounding disc as shown in Fig.l. The radii ri n 
and r ou t are the radii of inner and outer boundaries of the 
MC region, respectively, and Q\ and 62 are the corresponding 
poloidal angular coordinates on the horizon. 

(ii) The toroidal profile of the magnetic field is non- 
axisymmetric on the BH horizon, being expressed by a func- 
tion / ((/>) °f the azimuthal coordinate (j> as shown in Fig. 2. 
The magnetic field is assumed to vary in a power-law as 
the radial coordinate of the disc as suggested by Blandford 
(1976). 

(iii) The disc is both stable and perfectly conducting, 
and the closed magnetic field lines are frozen in the disc. 

(iv) The disc is thin and Keplerian, locating in the equa- 
torial plane of the BH with the inner boundary at the last 
stable circular orbit. 

Generally speaking, the magnetic field on the horizon 
is non-axisymmetric. Unfortunately, a reliable profile of the 
magnetic field on the horizon has not been given due to 
lack of knowledge about it. As a simple analysis we assume 
the profile function / (cj>) as shown in Fig. 2(a), which is ex- 
pressed by 



fW = 



1 + 5, < cf> < A(f>; 



(1) 



1, 



A(f> sC 4> sC 2tt. 



where y/ (Bjj) is root-mean-square of the magnetic field over 
the poloidal angular coordinate from 9\ to 82- As shown in 
Fig. 2(a), the parameter 8 is used to describe the strength of 
the bulging magnetic field in the region < <f> < A<f>, which 
is the azimuthal angular range of the hot spot on the disc. 
The drawback of the pattern given in Fig. 2(a) lies in the 
uncontinuity of the magnetic field at <j> = and Acj>. In fact, 
function / (</>) can be modified by a smooth joint-function as 
shown in Fig. 2(b), and we shall modify the model of rotating 
hot spot by using a more realistic profile function in future. 

The radial range of the hot spot is determined by a 
mapping relation between the angular coordinate on the BH 
horizon and the radial coordinate on the disc as follows: 



G(a»;^,n)d^, 



where 
G(a„;£,n) 



5 1 - 



(2) 



(3) 



2^(l + a ; x ^+2a5x™ 6 a )(l-2x^?- 1 +^x^e- 2 ) 



The angle 62 — 7r/2 is assumed in derivation. Equation 
associated with equation © is the revised version of the 
mapping relation derived in WXL. Parameter £ = r/r ms is 
a dimensionless radius of the disc in terms of the radius of 
the last stable circular orbit, r ms — M^ 2 ma , and a* = J / M 2 
is dimensionless parameter of the BH spin defined by the 
BH mass M and angular momentum J. Following Blandford 
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Figure 3. The curves of £ c (dot-dashed line) and £out versus a* for < a, < 1 with n = 1.1, 1.5 and 3 in solid, dashed and dotted lines, 
respectively. The values of f?i are given as (a)6*i = 7r/12, (b)f?i = 7r/6 and (c)$i = n/4. 



(1976) the normal components of magnetic field, B z , on the 
disc is assumed to vary with £ as follows. 



(4) 



In deriving equations 10 and (0 the boundary condition is 
used as given in WLM, i.e. 



2-Kr H B± = 2-kvju (r ms ) B z (r ms ) , 



(5) 



where B± is the normal component of magnetic field at hori- 
zon, and vjd {r ms ) is the cylindrical radius at r m3 , 



kud (r ms ) = MxLs V 1 + Xmta* + ^Xmsdl 



(6) 



Equation (0 implies the conservation of magnetic flux 
corresponding to the two loops of the same width, and the 
ratio of B± to B z (r ms ) varies from 1.8 to 3 for < a* < 1. 
Considering that the strength of the magnetic field at the 
horizon is likely greater than that in the disc by numerical 
simulation (Ghosh & Abramowicz 1997 and the references 
therein), we think that the relation (0 is more reasonable 
than the boundary condition B±_ = B z (r ms ) given in WXL. 

Substituting 9 = 9\ and £ = £ out into equation (0, 
we obtain an integral equation for determining the outer 
boundary £ ou t = r ou t/r ms of the MC region as follows: 



cos til 



G(a.;£,n)d$. 



(7) 



From equation we find that £ „t (o*;n, 9i) behaves as a 
non-monotonic function of a, for different values of n and 
8i as shown in Fig. 3. 

In WLM the transfer direction of energy and angular 
momentum between a rotating BH and its surrounding disc 
is discussed in detail by using the co-rotation radius r c , 
which is defined as the radius where the angular velocity 
Q.D of the disc is equal to the BH angular velocity £Ih ■ The 
parameter £ c (a*) = r c /r ms is a function of a, determined 
by the following equation, 



where 



2(l + «) 



(8) 



Clo — 7 — r, 

M( x 3 + a.) ' 



2M(l+ 9 ) ' 



(9) 



In equation @ x = y r/M = y^Xms is a dimensionless 
radial parameter on the disc, and x ms is defined a-s \ms = 
yJr ms /M. The MC region is divided by r c into two parts: 
the inner MC region (henceforth IMCR) for 1 < £ < £ c and 
the outer MC region (henceforth OMCR) for £ c < £ < f ont . 
The parameter £ c decreasing monotonically with a* is shown 
by the dot-dashed line in Fig. 3, and the transfer direction 
of energy and angular momentum is described as follows. 

(i) For 0.3594 < a, < 1 the transfer direction is from 
the BH to the disc with £ c < 1. 

(ii) For ^ a* < a"" 1 the transfer direction is from the 
disc to the BH with £ c > £ otl t, where a° ut is the BH spin 
corresponding to the intersection of the curve £ c (a*) with 
the curve £ ou t (a», n, 6i) as shown in Fig. 3. 

(iii) For a° ut < a, < 0.3594 the transfer direction is 
from the BH to OMCR with £ c < £ < £ ut, while it is from 
IMCR to the BH with 1< £ < f c . 

Since a hot spot is produced by the energy transferred 
from the BH to the disc, its radial range is expressed by 



£ c < £ < tout, for a% ut < a, < 0.3594 (easel) 
K £ < tout, for 0.3594 < a* < 1 {easel I) . 



(10) 



From Fig.3 we find that the outer boundary £ ou t is very 
close to £i„ = 1 in both case I and case 77, provided that the 
power-law index n is small, such as n = 1.1, 1.5. This result 
seems insensitive to the values of 9i . Therefore the hot spot 
is confined to a small radial region very near to the inner 
edge of the disc. 



3 RADIATION ENERGY AND ROTATING 
FREQUENCY OF ROTATING HOT SPOT 

Based on conservation of energy and angular momentum the 
basic equations for BH evolution in the coexistence of disc 
accretion and the MC process are written as 
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dM/dt = E ms M D ~ Pmc, 
dJ/dt = L ms M D -T M c- 



(11) 
(12) 



where Pmc and Tmc are MC power and torque in the MC 
process, respectively. In WXL we derived Pmc and Tmc for 
axisymmetric magnetic field, by using an improved equiva- 
lent circuit as follows: 

ft 0(l-0)G(a,;t' ,n) 



= ±(B 2 H )M 3 a«(l + q)x 



f£ (l-/3)G(q.;e',n) , 
Jl 2 CSC 2 S-(l-q) a ? ■ 



(13) 



(14) 



Equations 1131 and 1141 can be extended to the case for 
non-axisymmetric magnetic field as follows: 



P:naxis j^axis 
MC — A-TMC ! 

Tinaxis \rncixis 
MC — A+MC i 



(15) 
(16) 



where the superscripts 1 axis' and 'naxis' indicate the quanti- 
ties for axisymmetric and non-axisymmetric magnetic field, 
respectively. The parameter A is expressed by 



\ = {(l+S)e + (l-e)} 2 = (l + fe) 2 



(17) 



where e = A<j)/2n is a parameter indicating the azimuthal 
region of the bulging magnetic field on the horizon expressed 
by equation Q. Equations 1151 and 1161 can be derived in 
the same way as in WXL, provided that the magnetic flux 
between two adjacent magnetic surfaces A\& = Bh2tt-ouAI 
is replaced by 



A* = y/{Bjj)zu [(1 + S)A<j,+ (2tt - A4>)] Al 



(18) 



yf(Bl)2TTzu (1 + 5e) 2 Al, 



where equation (1181 is given for the bulging magnetic field 
expressed by equation Q. Obviously, equations 1151 and 
1161 reduce to axisymmetric ones for either 8 — ► or e — > 
with A = 1, and they also become axisymmetric ones for 
e -t 1 with A = (1 + 5)' 2 . Based on equations (1131 — 1161 
we can express the rate of extracting energy and angular 
momentum from the BH to the hot spot as follows: 



Pmc — (A 



Tmc 



1) Pmc + sPmc = (A + e - 1) P MC , 



(A 



1 \ r-i-iCLXlS 

i-J-l-MC ■ 



rjifixis 

SJ-MC 



(A + e- 



1 \ rriaxis 

1 ) 1 MC ■ 



(19) 
(20) 



From equations (1151 . 1161 . 1191 and 1201 we find that the 
ratio of P§$ to P^/a " is the same as that of T^ to Tl 
i.e. 



rnaxts 
MC 



R (5, e) 



nHS 


_ T MC 


A+e-1 


P MC 


5 T »i <■ . / ., 
MC 


A 


= 1 - 


1-e 





(21) 



where the ratio R (5, e) is used to express the fractional en- 
ergy and angular momentum attributable to the hot spot 
in the MC process. The curves of R (5, e) varying with the 
parameters 8 and e are shown in Fig. 4. 

From Fig. 4 we find that the ratio R (S, e) is indepen- 
dent of the parameters M, a* and n, and it increases almost 
linearly with the parameters 8 and e. For an appropriate 
angular region of the bulging magnetic field on the horizon, 
such as e — 0.1, we have the ratio varying from 0.1 to 0.18 
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Figure 4. The curves of R(8, e) (a) versus 8 for < 8 < 0.5 with 
e = 0.1,0.05 and (b) versus e for < e < 0.1 with <5 = 0.5,0.3 in 
solid and dashed lines, respectively. 



for a small value range of the strength of the bulging mag- 
netic field, < 8 < 0.5. These results show that the effects 
of the bulging magnetic field on the energy and angular mo- 
mentum transferred to the hot spot are significant. 

Since the hot spot is located very near to the inner edge 
of the disc, it should be optically thick. So we can define 
its effective radiation temperature (Tijs) e ff as a blackbody 
spectrum: 



>eff 



[(Fda- 



F M S c)/c 



,1/4 



(22) 



where a is the Stefan-Boltzmann constant. The fluxes Fda 
and Fmc are radiated from the hot spot due to disc accre- 
tion and the MC process, respectively, and are expressed as 
follows (Page & Thorne 1974; Li02a): 



F DJ 



F 1 



f (E+-Q D L+y 2 x 

(E+ -n D L + )HHs£d£. 



dSl 



(23) 



(24) 
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In this paper we only discuss the hot spot produced by 
the bulging magnetic field on the BH with spin greater than 
a» = 0.3594, which corresponds to case II in equation 1101 . 
In equations 1231 and 1241 E + and L + are specific energy 
and angular momentum of accreting particles, respectively, 
and are expressed respectively by (Novikov & Thorne 1973) 



E+ = 



(l-3x- 2 + 2a*x" 3 ) 1/2 
and 

r+ M X (l-2a, X - 3 



(25) 



(26) 



(l-3 X - 2 + 2a* X - 3 ) 1/2 ' 

In equation 1241 1 Hhs is the flux of angular momen- 
tum transferred from the BH to the hot spot by the bulging 
magnetic field, being related to T^q by 

A-ketHhs = 2rA</>H HS = dT^c / dr (27) 

Combining equations ©, 1141 . 121 It and 1271 . we have 



Hhs/H 



A(a„£,n,5,e)£ n , 1 < £ < £ oui 

o, £ > £ out 



(28) 



where 



. , c c \ (l+2S+8 2 e)a,(l-0)(l+ q ) 

A(a*,$,n,5,s) = \ g _ . F A {a.,0 



1 F A (a*,0 



(29) 



^/(l + ai X - 4 B +2* 2 ,X^ s )(l-2x;. 



and 



Ho = (B 2 H ) M = 1.48 x 10 13 x B 2 m B H g ■ s~ 



(30) 



In equation 1301 1 B4 and tubh are yj (Bjj) and M in the 
units of W 4 Gauss and one solar mass, respectively. 

Since the magnetic field on the horizon is brought 
and held by its surrounding magnetized disc (Macdonald 
& Thorne 1982; Thorne, Press & Macdonald 1986), there 
must exist some relations between Bh and Md- As a mat- 
ter of fact these relations might be rather complicated, and 
would be very different in different situations. One of them 
is given by considering the balance between the pressure of 



the magnetic field on the horizon and the ram pressure of 
the innermost parts of an accretion flow (Moderski, Sikora 
& Lasota 1997), i.e. 

B 2 „/{8n) = P ram ~ pc 2 ~ M D /(4nr 2 H ), (31) 

From equation 1311 we assume the relation as 

M D = (B 2 H ) M 2 (1 + g) 2 /2. (32) 

The accretion rate Md is still axisymmetric in spite of 
the existence of the bulging magnetic field on the horizon, 
provided that the extra angular momentum transferred into 
the disc is totally radiated away from the disc. Incorporating 
equations lESl . t2H and 1281 . we have 

Fda/Fq = /da (a*,£) 

/^(a*,0=- i ^ [ J±g^ZTf X (33) 
and 

F mc/ F = /MO, 

}mc 



MdClp 



(E+-n D L+) 2 ' « d « " 



(34) 



with 



F = (B 2 H ) c = 2.998 x 10 18 B 2 erg ■ cm' 2 ■ s" 1 . 



(35) 



Inspecting equations 1331 — (13511 . we obtain the follow- 
ing results: 

(i) Both Fda and Fjgg are proportional to (B^), and 
are independent of the BH mass due to equation 1321 . 

(ii) The bulging magnetic field on the horizon only af- 
fects Fmc rather than Fda, since Md is independent of the 
bulging magnetic field. 

By using equations 1331 and 1341 we have the curves of 
Fda and F^c varying with £ as shown in Fig. 5. 
From Fig. 5 we obtain the following results: 
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Table 1. Radiation energy and frequency of a rotating hot spot produced by the MC process and disc accretion with B4 = 1, 8 = 0.5 
and e = 0.2. 



0* 


n 


( T sr) eff ( K ) 


£g;| x (kev) 


(max 






1.1 


7.722 xlO 4 


6.661xl0~ 3 


1.348 


2620.37 


0.6 


1.5 


7.666 xlO 4 


6.612xl0~ :i 


1.367 


2569.23 




3.0 


7.440 xlO 4 


6.417xl0~ 3 


1.447 


2366.91 




1.1 


1.038 x10 s 


8.953X10" 3 


1.380 


3657.66 


0.8 


1.5 


1.029 x10 s 


8.875xl0- :i 


1.396 


3599.78 




3.0 


9.950 xlO 4 


8.582x10^3 


1.438 


3455.97 




1.1 


1.994x10 s 


1.720X10" 2 


1.250 


11061.7 


0.998 


1.5 


1.982x10 s 


1.709X10- 2 


1.247 


11084.3 




3.0 


1.937x10 s 


1.671xl0~ 2 


1.227 


11261.7 



(i) The relation F^c > Fda always holds for 0.3594 < 
a* < 1, implying that the contribution from the MC process 
always dominates over that from disc accretion for the above 
value range of the BH spin. 

(ii) Both F^q and Fda approach zero as £ — » £ in = 
1, and both vary non-monotonically with £. The flux F^q 
varies much more sharply than Fda does, attaining its peak 
value between the inner and outer boundaries of the hot 
spot. 

(iii) The peak of Fj^c ls closer to the inner edge of the 
disc than that of Fda • The more is the value of a*, the 
greater is the peak value, and the closer is the peak to the 
inner edge. 

Substituting equations 1331 and 134H into equation 1221 . 
we have the effective temperature and the corresponding 
radiation energy of the hot spot as follows: 

(T H s) eff = T [f D A + /mc] I/4 , (36) 

E HS = k B (T HS ) eff = E Q [f DA + f M c] 1/4 , (37) 
where 

To = {F /a) 1/4 « 4.8 x 10 5 B 4 1/2 A", (38) 
and 

E = k B (F /a) 1/4 « 4.14 X 10" 2 B\ /2 kev. (39) 

The maxima of (Tns) a ff and Ehs corresponding to the 
different values of the BH spin and the power-law index are 
shown in Table 1. 

It is easy to prove that the maxima of either (Tffs) 
or Ehs increase monotonically with the increasing BH spin. 
For the fast spinning BH with spin 0.6 < a* < 0.998 
we conclude that the radiation energy of the hot spot 
reaches 0.006 ~ Q.QUkev, provided that the bulging mag- 
netic field reaches B4 or 10 4 Gauss, while this value range 
becomes 0.6 ~ 1.7kev for the bulging magnetic field reach- 
ing 10 4 B4 or W a Gauss. So the hot spot can have energy 
reaching soft to hard X-ray in the MC process, provided 
that a fast-spinning BH has the bulging magnetic field of 
10 4 ~ 10 s Gauss on the horizon. 

Since the magnetic field lines are frozen in the disc 
plasma, the rotating frequency of the hot spot is regarded 
as QPO frequency. As suggested by Nowak and Lehr (1998), 
the QPO frequency can be worked out by calculating Ke- 
plerian angular velocity at the place where r 2 F attains its 




0.4 0.6 0.8 1 

a* 



Figure 6. The contours of Ehs of constant values (kev) for 
0.3594 < a* < 0.998 and 10 3 < B 4 < 10 4 with 6 = 0.5, e = 0.1 
and n = 1.1,1.5, 3.0 in solid, dashed and dotted lines, respec- 
tively. 



maximum, and F is the sum of Fmc and Fda- We can define 
a dimensionless function as follows: 

Fqpo = r 2 F/r 2 ma F = £ 2 (Fba/F + F^c/Fo) (40) 

Incorporating equations 13311 . 13411 and 140H . we obtain £ max 
corresponding to the maximum of function Fqpo, and the 
QPO frequency vhs can be calculated by 

vhs = ^KmixL + a*)~\ (41) 

where we have vhs = S7u/27r| ?=? and vq = (msff) _1 X 
3.23 x 10 4 Hz. The values of £ max and vhs corresponding to 
different values of BH spin and power-law index are shown 
in Table 1, where £ max turns out to be very close to the 
inner edge of the disc. For the fast-spinning BH with spin 
0.6 < a» < 0.998 we have vhs varying from ~ 2AkHz/m,BH 
to ~ HkHz/nriBH with tubh > 3. 

In order to describe the variation of Ehs and vhs of the 
rotating hot spot more clearly we plot the contours of Ehs 
and vhs in terms of the concerning parameters as shown in 
Figs 6 and 7, respectively. 

From Fig. 6 we find that less magnetic field is required 
as the BH spin increases for the given value of Ehs- Further- 
more the more is the BH spin, the more steeply the contours 



An analytic model of rotating hot spot and kHz QPOs in X-ray binaries 7 



30 



25 



20 



15 



10 






0.02 




0.01 







> 










- 0.01 












- 0.02 




- 0.03 




- 0.04 



0.6 



0.4 



0.6 



0.8 



0.7 0.S 

a. 



0.9 



Figure 7. The contours of vhs of constant values (Hz) for 
0.3594 < a» < 0.998 and 3 < m BH < 30 with <5 = 0.5, e = 0.1 
and n = 1.1,1.5, 3.0 in solid, dashed and dotted lines, respec- 
tively. 



of Ehs decline. This result implies that the MC process be- 
comes very efficient in transferring energy from a very fast 
spinning BH to the hot spot. 

From Fig. 7 we find that more BH mass is required 
as the BH spin increases for the given values of vhs- Fur- 
thermore the more is the BH spin, the more contours of 
vhs will be crossed over for the given BH mass. This result 
implies that vhs increases very rapidly as the BH spin o* 
approaches unity for the given BH mass. 

The BH spin acts as a bridge to link the contours of 
Ehs and vhs in Figs 6 and 7. This linkage is helpful to fit 
observations of kHz QPOs in X-ray binaries in some scope. 



4 FLUCTUATION OF PARAMETERS AND 
WIDTHS OF QPO FREQUENCIES 

ft is shown that the widths of observed frequencies of kilo- 
hertz QPOs in X-ray binaries are rather narrow, usually 
varying from several to some dozens of hertz, e.g. as shown 
in Tables 1-3 given by van der Klis (2000). However, it is 
difficult for us to explain the narrow range of QPOs by the 
rotating hot spot produced by the stationary bulging mag- 
netic field on the horizon. 

As a matter of fact we can fit the narrow range of QPO 
widths by the fluctuation of S and e in our model. The ar- 
gument is given as follows. From equation 1411 we find that 
QPO frequency is determined by £ ma x, and the latter de- 
pends on the maximum of function Fqpo in equation 1401 . 
which consists of the sum of F D a/Fo and F$%/F . As men- 
tioned above that the peak of F^c / Fa is located closer to 
the inner edge of the disc than that of Fda/Fo, so we infer 
that the narrow range of QPO widths might be produced 
by the variation of Fmc/Fo, provided that the parameters 

5 and e are fluctuating rather than constant. The curves of 
dvns I vhs versus a* for the given fluctuation of 8 and e are 
shown in Fig. 8. 

From Fig.8 we have the following results, 
(i) The value of dvns/vns varies between -0.04 to 0.02 
with |A<5| = 0.3 and Ae = for 0.54 < o» < 0.998. This 



(a) 



0.001 



0.0005 



0.0005 



0.001 



0.6 



0.7 



0.8 



0.9 



(b) 

Figure 8. The curves of cLvhs/vhs versus a* for 0.54 < a* < 
0.998, (a) n = 3, e = 0.2, Ae = 0, 5 = 05, AS = 0.3 (solid line), 
AS = -0.3 (dashed line), (b) n = 3, 5 = 0.5, AS = 0, e = 0.2, 
Ae = 0.1 (solid line), Ae = -0.1 (dashed line). 



result implies that the range of QPO widths is of dozens 
of hertz, and could be produced by the fluctuation of the 
strength of the bulging magnetic field of the BH as shown 
in Fig.8(a). 

(ii) The effect of the variation of e on the range of kHz 
QPO width is about one order of magnitude less than that 
of 8 as shown in Fig. 8(b). 

The values of dvns / vhs corresponding to different val- 
ues of the concerning parameters are listed in Table 2. As 
shown in Fig. 8(a) and Table 2 we can fit the range of kHz 
QPO width by a variety of parameters in our model, such 
as a*, n, e, 8, Ae and AS. 



5 SUMMARY 

In this paper an analytic model of rotating hot spot is pro- 
posed to explain kHz QPOs in X-ray binaries by using the 
MC of a rotating BH with the disc. Six parameters are used 
in our model. 

(i) The parameters S, e and B4 are given for describing 
the bulging magnetic field on the horizon by equation Q. 
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Table 2. The values of dvus/^HS corresponding to different values of the concerning parameters with e = 0.2, Ae = and S = 0.5. 



a* 


n 








dvnsl v HS 


















AS = 0.1 


AS = -0.1 


AS = 0.2 


AS = -0.2 




1.1 


2620.37 


2.42X10" 3 




4.46X10" 3 




0.600 


1.5 


2569.23 


3.01X10" 3 


-3.63X10" 3 


5.55X10" 3 


-6.36X10" 3 




3.0 


2366.91 


7.35x10-3 


-8.93xl0- 3 


13.5X10- 3 


-20.0X10" 3 




1.1 


3657.66 


1.94xl0- 3 


-2.34XKT 3 


3.57X10" 3 


-5.23X10" 3 


0.800 


1.5 


3599.78 


2.41xl0~ 3 


-2.91xl0~ 3 


4.43xl0~ 3 


-6.50xl0~ 3 




3.0 


3455.97 


5.07X10" 3 


-6.17X10" 3 


9.32X10" 3 


-13.8X10" 3 




1.1 


11061.7 


1.31X10- 3 


-1.58XKT 3 


2.42X10" 3 


-3.53X10" 3 


0.998 


1.5 


11084.3 


1.44xl0~ 3 


-1.74xl0~ 3 


2.65X10- 3 


-3.87xl0~ 3 




3.0 


11261.7 


1.78xl0~ 3 


-2.15X10" 3 


3.27X10" 3 


-4.81X10" 3 



Specifically speaking, the parameter B4 is used to indicate 
root-mean-square of the magnetic field over the poloidal an- 
gular coordinate from Qi to #2. The parameters S and e are 
in charge of the fractional energy and angular momentum 
attributable to the hot spot by equation 1211 1. The fluctua- 
tion of S and e can affect the range of kHz QPO widths. 

(ii) The parameters ttlbh and a* are given for describing 
the mass and spin of the BH, which determine the frequency 
of QPOs by equations 1411 , 

(iii) The parameter n is used to indicate the concentra- 
tion of magnetic field on the central region of the disc, and 
it is involved in equation Q to confine the radial width of 
the hot spot. In addition, it is involved in equation 1401 to 
determine £ max . We can find that the index n plays a role 
of fine-tuning energy and frequency of the hot spot and 
the range of kHz QPOs widths as shown in Table 1, Fig. 6, 
Fig. 7 and Table 2, respectively. 
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